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Abstract. The Hamiltonian representation of magnetic field lines enables to study 
their confinement properties in tokamaks through the use of symplectic maps such 
as the symmetric tokamap and its bounded version, the latter being introduced 
here. In this time-independent purely magnetic framework, we observed the drastic 
improvement on the confinement of magnetic field lines produced by the local vanishing 
of the shear profile. This amounts to a non-twist condition that notably acts in the 
same way the safety profile being (non-strictly) monotonic or having reversed-shear. 
We single out the effect of the amount of flatness of the safety profile in the vicinity of 
its zero shear point. All other things being equal, the beneficial effect of the vanishing 
of the shear profile is shown to be increased if the radial extent of the low-shear region 
is increased. To be specific, the low-shear region induces the formation of a belt of 
robust KAM tori acting as an internal transport barrier whose width is all the larger 
as the extent of the low-shear region is broad. 
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1. Introduction 

In order to approach ignition, it is required that the present hot tokamak plasmas have 
a better energy confinement. The energy confinement time is commonly limited by 
deleterious small-scale turbulent transport. Because internal transport barriers (ITB) 
[H El [3] reduce or even quench the turbulent transport within the inner plasma region 
improving energy confinement in the plasma core which should yield high fusion gains, 
they have become an essential ingredient of advanced scenarios in which the tokamak 
could operate as an economic steady-state fusion reactor. 

The physics of the triggering and sustainment of ITB as well as the explanation 
of their very nature are still subjects of intensive investigations [3]. Besides an 
increasingly detailed empirical and experimental knowledge, their theoretical modeling 
is still incomplete. It is commonly admitted that two ingredients should control their 
formation, namely the E x B fiow shear and the magnetic shear through the safety factor 
(g) profile. An increased number of works have recently reported the importance of the 
g-profile. In particular, Eriksson et al. [1] demonstrated experimentally in the JET 
tokamak the critical role of the safety factor profile in the formation of ITB by using 
it as a control parameter in discharges where other parameters, including the E x B 
fiows, were kept the same. They obtained an ITB for a non-monotonic q profile (i.e. 
with reversed shear) but not with a monotone q. Actually these reversed shear scenarios 
have been recognized as favorable experimental conditions to obtain ITB whereas the 
ITB's locations have been mostly related to the low shear zones. 

In this article, we wish to study the appearance and robustness of ITB within a 
purely magnetic and static (time-independent) framework following a line of approach 
notably initiated by Balescu and coworkers. This simplified, yet relevant, approach is 
guided by the fact that, in a first approximation, charged particles in a tokamak follow 
the magnetic field lines so that transport properties for the laters transfer to the formers. 
We shall focus on the infiuence of the behaviour of the g-profile in the vicinity of its 
vanishing shear region on the appearance and robustness of ITB. To be more specific, we 
wish to single out the role and importance of the amount of flatness of the g-profile in 
the vicinity of its minimum (for the reversed shear case) or in the vicinity of its inflexion 
point (in the case of a non-strictly monotonic profile). 

In Section [21 the Hamiltonian formulation of magnetic field lines is recalled and 
symplectic maps are introduced. In order to capture more accurately the magnetic 
field lines, we shall implement the symmetric tokamap and its bounded version, instead 
of the standard tokamap. In Section [3], we shall study numerically the infiuence of 
the flatness of the g-proflle, in the vicinity of a given inflexion or minimum point, on 
the overall magnetic fleld lines stochasticity. To do so, we shall consider a family of 
polynomial winding number proflles. In Section HI we eventually consider a sample of 
tokamak-oriented realistic g-proflles with reversed shear. This brings another evidence 
on the favorable role played by an increased extension of the low-shear region around 
the minimum of the g-proflle. We conclude in Section [51 
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2. Hamiltonian formulation for magnetic field lines and maps 

2.1. Hamiltonian framework 

In order to investigate the role of the g-profile on the existence and robustness of ITBs, 
we shall use the Hamiltonian structure of the magnetic field lines. We shall use in 
the following traditional generalized toroidal coordinates ip, 6, ( where ip is the flux 
coordinate (so that ip = and = 1 correspond to the magnetic axis and to the the 
plasma boundary, respectively), 9 is the poloidal angle and ( the toroidal angle. From 
Maxwell's equation V.B = 0, the "equations of motion" for magnetic field lines j5] in a 
tokamak are 

dip _ _dF dO _dF 

~dC ~ ~~de ' dC~d^ ^ ^ 

where F{iIj, 6, () is the poloidal fiux. Equations ([1]) have a Hamiltonian structure: 
the toroidal angle ( plays the role of "time", F is the Hamiltonian and 6 and ip are 
canonically conjugated coordinates. In order to consider perturbations from the ideal 
i.e. integrable case, where F = -Fo(^) and magnetic surfaces are nested tori wound 
around the magnetic axis, let us introduce a stochasticity parameter, K, such that the 
Hamiltonian F is given by the sum of an unperturbed term Fq and a perturbation 
6F{^,e,0 : 

F = Fo{^) + K5Fi^,eX). (2) 

As K increases, magnetic surfaces may no longer exist and the problem of their existence 
is formally equivalent to that of the integrability of a Hamiltonian with one-and-a-half 
degrees of freedom. In realistic tokamaks, these magnetic perturbations may result 
e.g. from coil imperfections and/or internal instabilities (tearing). Within the above 
formalism, the safety factor q{ip) is simply related to the unperturbed Hamiltonian since 
qlip) = l/Wltp) where W{iIj) = dFo/dip is the winding number. Therefore q{ip), the 
ratio of the number of toroidal turns per poloidal turn, is the inverse of the unperturbed 
Hamiltonian. Moreover, s = dlnq/d\nip is the shear profile. 

2.2. Maps 

2.2.1. Tokamap Instead of performing a time consuming integration of the field line 
differential equations ([I]), we shall use a simplified but relevant approach, namely the 
Poincare map associated to a given poloidal cross-section [6]. It is defined from the 
intersection points of a magnetic field line starting at position [i/jq, 9q) with some given 
poloidal section C, = est. The intersection point after i/ toroidal turns is denoted by 
{ipu,(^u)- The map we shall be interested in, the tokamap, was introduced by Balescu 
and coworkers [Zl El [9] . Its implicit form is 

ip^+i =ipu- K ^"^^ sin 9^ 

1 + Vu+l 

(3) 
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K 1 

Explicitly, we have 

= I (^('«'- + VlPU'v.o.W + H'v) (4) 



where 



P{iP,e) = ip-l-Ksme. (6) 

This map is an area-preserving map which is compatible with the toroidal geometry [7] 
due to the following important physical properties : (i) the toroidal flux ip is always 
positive (this is necessary since ■?/' ~ where r is the minor radius of the tokamak) and 
(ii) the magnetic axis is invariant [tp = is mapped to = 0). It has proved to be 
an important model since it correctly describes the qualitative features of the magnetic 
field lines stochasticity known in tokamak physics. 

2.2.2. Symmetric tokamap Some studies have recently pointed some limitations of 
this model. It actually turns out that despite a qualitative agreement with the actual 
magnetic field behaviour obtained from the continuous Hamiltonian system formed 
by Eqs. ([1]) and ([2]), the tokamap is far less accurate than its symmetric version as 
proposed by AbduUaev [IHIE]. We shall therefore implement in this study the symmetric 
tokamap instead of the standard one. The symmetric tokamap very closely describes 
the continuous Hamiltonian system ([2]) with 

^ ' ^ n=-M 

and 2M +1^1, from which it is derived. Its implicit form is given by 

K ^ 

e.+i = + 2'kW{^>,) (8) 

i^v+i = "^u+i - TTTTir, sm^^+i 

e.,, =e.,,-^^^^l-^cos^.,, 

where K is the stochasticity parameter and the auxiliary variables and ^ are 
determined iteratively by the map from a given starting point {do, ipo)- Contrarily to the 
tokamap, a numerical approximation is necessary in order to solve the last equation of 
the symmetric tokamap with respect to variable 6u+i, but this is an affordable task 
with usual interpolation methods such as Newton's method or the Brent's algorithm 
[TT] that we used here. 
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2.2.3. Bounded symmetric tokamap Another adjustment to the tokamap has also been 
recently proposed: the bounded tokamap. Whereas the tokamap does not impose any 
upper boundary limit for ip, so that it can be directly used in the ergodic divertor 
problem [12] where magnetic field lines may go outside the last magnetic surface (defined 
as = 1), its bounded version is constructed by requiring that ip = 1 cannot be crossed. 
We used the approach developed in Refs. [131 E] to propose the symmetric version of 
the bounded tokamap with 

K 

^ ^ n=-M 

and 2M +1^1. This form of perturbation ensures that not only the magnetic axis, 
corresponding to ip = 0, but also the plasma boundary, corresponding to ip = ^ are 
invariant under the map and cannot be crossed. 

Then, applying the same method [lOl [6] as in Sec. 12.2.21 we propose the symmetric 
version of the bounded tokamak as 







K 

ZTT 

^(1 -2^,)cos^, 

ZTT 




= 0,- 


e.+i 


= e,+ 


27rW^(^^) 






K 

- — ^^(l-^^)sin^^+i 

ZTT 










K 

- — (1 -2^^)cos^^+i. 
In 



The first equation of flTUl) has two roots for \E'j, and the positive one is retained. This 
modeling is valid for not too large values of the stochasticity parameter as it is required 
that K/271 < 1. 



2.3. Motivations 

From the theoretical point of view, the effect of various parameters describing the q- 
profile, such as the boundary values of the g-profile, the value of the minimum of the q- 
profile in the reversed shear version (in the so-called rev-tokamap), for Hamiltonian maps 
have already been considered [T^l US]. Yet, this has been mainly limited to restricted 
choices of q profiles. In addition, to the authors' knowledge, an essential point has 
not been investigated systematically: the effect of the local flatness of the g-proflle on 
(hopefully) realistic models for magnetic fleld lines in a tokamak. This is in our sense 
a crucial point since there is strong experimental evidence that the radial extent of the 
low shear region of safety factor's profile plays an outstanding role on ITBs' appearance 
and robustness [I71[2],[3]. Thus, modeling is needed in order to asses one of the ''critical 
physics issues relevant to the extrapolation of ITB regimes to next-step experiments, 
such as ITER " [18] : the optimum g-profile. This forms the ultimate perspective of 
the present work. 
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3. Numerical results on the stochasticity of magnetic field lines for 
polynomial winding number profiles 

We are interested in assessing the influence of the extent of the low shear region on 
magnetic conflnement. Since it happens that the Hamiltonian ([2]) is more directly linked 
to the winding number, since W{ip) = dFo/dip, than the g-proflle, we shall consider it 
in the following. We are then particularly interested in the influence of the behaviour 
of the PF-proflle in the vicinity of some unique ipo such that W'lipo) = 0. If such a ipo 
exists, the twist condition for the corresponding map is violated, that is the symmetric 
tokamap is non-twist. A straightforward way to study such behaviour is to consider the 
family of polynomial PF-profiles : 



where n is a positive integer representing the degree of degeneracy of the zero ipo of the 
shear profllelj associated to ffTTl) . Although there is an extensive bunch of experimental 
evidence that the maximum value of the winding number influences 1TB formation in 
the reversed shear case p!7l El [3] , this will not be our point here and we shall rather 
impose here on an arbitrary basis Wq = 5/2it, c = 1.3 and ipQ = 0.5 in order to focus on 
the flatness of the W^-proflle only. 

For every proflle (|TT|) . with the degree n going from 1 to 8, the phase portrait of 
the symmetric tokamap (IHl) has been computed. The results of this iterative process, for 
ip in the physical domain [0, 1], are given for the same value = 3 of the stochasticity 
parameter in Fig. [H Several comments are called for. First of all, it is clear from 
the flgure that the case n = 1 is dramatically different from the others. In the case 
n = 1, the map is twist as W'lip) does not vanish. The flgure exemplifles the dramatic 
improvement on confinement induced by the non-twist condition that is satisfled when 
n > 2. This has been discussed recently on a general basis by Rypina et al. 
Secondly, when the non-twist condition is satisfled, the core magnetic fleld lines are 
almost integrable. The central non-chaotic region is all the broader in ^/^-space as the 
ly-proflle is flat around ipQ, namely as n is large. Yet this effect saturates rapidly: for 
n > 6, the phase space portraits are almost indistinguishable. This is not surprising 
as the corresponding g-proflles are indeed almost identical. Thirdly, there is a perfect 
progression in the phase space appearances between the odd and even n-exponents as n 
increases in the non-twist case i.e. for n>2. This means that the fact that the g-proflle 
has a reversed shear or not is not important in terms of the conflnement of the magnetic 
fleld lines. What matters is the degree of degeneracy of the zero of the shear proflle or, 
so to speak, the amount of local flatness of the l^-proflle. To check this explicitly, let 
us consider the two following proflles 



W{%Ij) 



(11) 



c{iP - V^o)' + Wo 

■c\{^-^o)? + Wo 



(12) 



(13) 



\ More precisely, -00 is a zero of the shear profile when n > 2. 
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Figure 1. Phase space portraits of the symmetric tokamap ([8]) for different values 
of the order of degeneracy n of the profile (fTTj) . Parameters are K — 3, Wq = 5/2tt, 
c = 1.3 and = 0.5. 

where Wm f|T2|) and Wrs ( IT3|) correspond respectively to a monotonic and to a reversed- 
shear profiles with the same local flatness (n = 3). The results are displayed in flgureO 
This shows that the global behaviour of the magnetic fleld lines is indeed the same in 
both cases: the outer stochastic region is identical. There are only some limited apparent 
topological discrepancies in the inner ordered zone. This third point is rather novel since 
previous studies [15] commonly used to focus solely on reversed- shear g-proflles to get 
the non-twist condition. 

We have brought some numerical evidence of the influence on magnetic conflnement 
of the amount of flatness of the g-proflle around its extremum. In order to deflnitely 
ascertain this phenomenon, we eventually wish to avoid the possible bias induced by 
using different edge values of W for each n in Fig. [1] It actually happens that increasing 
the edge value of a polynomial VT-proflle ( ITTi) with given degree n may induce the 
stochastization of the, otherwise regular, core magnetic field lines. This is illustrated in 
Fig. [3] for the case n = 2. Therefore, in order to single out the effect of the flatness of 
the IV-proflle i.e. of the radial extent of the low-shear region around an extremum, the 
study shown in Fig. [1] was revised for the polynomial proflles flTTj) with variable c = c„ 
such that W{0) = W{1) be a given constant for all n. Results are depicted in Fig. Hlfor 
even values of n between 2 and 8. This flgure shows that, all other things being equal, 
the radial extent of the regular magnetic field lines is all the wider as the low shear 
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0.9 - ■ - ■ , ' 0.9 - 




Figure 2. Phase space portraits of the symmetric tokamap jS]). (left) Wm{ip) is given 
by dni); (right) WVs(V') is given by Parameters K = 3, Wq = 5/2tt, c = 1.3 and 

V-o = 0.5. 



region around the extremum of the g-profile is broad. 

Finally, we wish to insist on the fact that the above conclusion does not depend on 
the particular choice of the map. It remains valid when considering the - less accurate - 
standard tokamap and the bounded version of the symmetric tokamap. In this case, the 
appearance of internal transport barriers can be seen more easily. This can be shown in 
Fig. [5l Whereas global chaos can be observed for n = 1, the non-twist effect that takes 
places for n > 2 dramatically improves confinement. A clear internal transport barrier 
emerges for n = 2 separating the internal and the outer chaotic regions. Increasing n, 
that is the extent of the low shear region, induces a drastic reduction of chaos. 

4. Influence of the radial extent of the low-shear region for some realistic 
reversed shear g-profiles 

The previous academic case of polynomial centered VT-profiles is not realistic. We wish 
in this Section to briefly test our results on the beneficial effect of the increased radial 
extent of the low shear domain on more experimentally compatible g-profiles. We shall 
here consider reversed shear profiles leading to ITBs. Typical locations of the minimum 
of the g-profile correspond to a normalized radius of 0.3-0.4 [1] and to a central g(0) ~ 3 
while qmin, the minimum value of g, is about 2. We have considered three g-profiles 
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Figure 3. Phase space portraits of the symmetric tokamap ([8]) with W{ip) = 
5/27r + c{4> - .5)2 and X = 3 for: (left) c = 2, (right) c = 1. 



compatible with these reahstic data and computed the corresponding phase portraits 
for the magnetic field lines. The results are given in Fig. [61 It is once again clear from 
this figure that, for increasing fiatness of q about its minimum and for constant K, the 
extent of the non-chaotic, regular, magnetic field lines in the plasma core is broader. In 
order to show that this effect is robust, namely map independent, we have computed the 
phase space portraits for the symmetric bounded tokamap ( fTOj) for the same g-profiles. 
Results are plotted in Fig. [61 We do observe once again that for increasing fiatness 
there is chaos reduction. 

A quantitative indicator of the degree of fiatness of some iy-profile around its ex- 
tremum located at ipo could be given by fwii'o) = limi/i^i/^o l^l"^) ~ W^ipo) \ / ^^li^ — 'ipo 
As shown previously on polynomial profiles, for a given stochasticity parameter, con- 
finement is improved as fwii'o) increases. However, fwii^o) gives the power of the first 
no n- vanishing term in the Taylor expansion of the ly-profile around tpo- Since realistic 
reversed-shear profiles are generically approximated locally by a parabola, this would 
yield /ly(V'o) = 2. This has been checked on the above profiles. It is thus necessary to 
produce a practical indicator of the local amount of low shear. Let us then assume that 
the reversed shear profile may be locally Taylor expanded as 

WW = WitPo) + ^W"MStP' + 0(5^3), (14) 
with Sip = ip — TpQ. Then, for a given Witpo), the local fiatness of W increases as 
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n=2;c=1 n=4;c=4 n=6;c=16 n=8;c=64 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



e/2ii 0/211 6/271 6/271 



Figure 4. (up) Reversed-shear q-profiles given by = l/[5/27r + Cn{ip — -5)"] with 
identical boundary values and increasing central flatness (with n even between 2 and 
8) with (down) their corresponding phase space portraits of the symmetric tokamap 
l^ior K ^ 3. 

|iy'(^/'o)| becomes smaller. This is in agreement with the results given in Fig. [6]where a 
diminution of W"{ipo) corresponding to an increased local flatness is observed from left 
to right. More precisely, one obtains numerically W"{ipo) = 10.6 for the profile located 
at the left of Fig. El W"{^o) = 7.8 for the middle profile and W^"(V'o) = 7.4 for the 
right one. Finally, let us note that in this Figure as well as in Fig. [5l the stochasticity 
parameter K, that is proportional to the relative amount of magnetic perturbation [8], 
has been taken unrealistically large. This points to the fact that the effect of low shear 
is much stronger than usual KAM results on the survival of tori that would be obtained 
for monotone g-profiles and corresponding twist maps [See e.g. the large scale chaos 
exhibited in the tokamap phase portrait of Fig. 9 in the work by Misguich et al. [8] 
that is obtained also for K = 6]. 

5. Conclusions 

The Hamiltonian formalism has been used to account for the behavior of magnetic 
field lines under some perturbations of given equilibrium safety factor profiles. The 
corresponding symplectic maps were the symmetric tokamap and its bounded version in 
which the last magnetic surface cannot be crossed. In this time-independent framework. 
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Figure 5. Phase space portraits of the bounded symmetric tokamap pO)l for different 
values of the order of degeneracy n (n = 1, 2, 3 and 4) of the profile (fTTj) . Parameters 
are K = 6,Wo = 5/2tt, c = 1.3 and ipo = 0.5. 




Figure 6. (top) g-profiles of increasing local flatness around Qmin', (middle) 
corresponding phase space portraits of the symmetric tokamap ([5]) for 500 iterations 
and K = 6; (bottom) same as above for the bounded symmetric tokamap PH)) . 
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we have shown the drastic improvement on the confinement of magnetic field fines 
produced by the local vanishing of the shear profile: the existence of some ipo such that 
g'('0o) = is a non-twist condition that acts in the same way the safety profile being 
monotonic or having reversed-shear. All other things being equal, the beneficial effect 
of this vanishing of the shear profile was shown to be increased if the radial extent 
of the low-shear region was increased. To be specific, the low-shear region induces 
the formation of a belt of robust KAM tori acting as an internal transport barrier 
whose width is all the larger as the extent of the low-shear region is broad. This 
conclusion has been derived here without any consideration of dynamical stability. It is 
important to note that there are additional dynamical arguments that favor a reduction 
of (electrostatic) turbulence for low magnetic shear, since it is well known that the linear 
growth rates of drift and ballooning type modes vanish at zero shear [20j . 

Finally, we wish to conclude by discussing the validity of this simplified purely 
magnetic approach. An obvious limitation of the model comes from its lack of self- 
consistency. In particular, there is no retroaction on the turbulence level. However, 
we would like to stress that such a simplified picture, being in a Hamiltonian form 
(due to the magnetic field being divergence- free), inherits from robustness properties 
w.r.t. fiuctuations. This is a general property of Hamiltonian systems and ensures that 
their long-time simulations are valid, at least in an average sense [21]. In the present 
model, this means that introducing some small perturbation in the symplectic symmetric 
tokamap, such a perturbation being e.g. a reduction of the magnetic turbulence [22] or 
another perturbation that may refiect the effect of electrostatic turbulence, the essential 
features of the model are fully preserved. 

It is now desirable to improve the connection of the present results with tokamak 
experiments and complete the modeling to address thermal transport. This would 
amount to couple the present magnetic approach to some heat transport modeling. 
Nevertheless, we believe that the benefits of the purely magnetic approach, initiated 
notably by Balescu, have not been exhausted yet. The role of the value of the q- 
profile in the low shear zone as well as the role of edge-gradients of the g-profile are 
other challenging important questions on which we shall be focusing our attention in a 
shortcoming future. 
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